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ABSTRACT 

We investigate the dynamical effects of rapid gas expulsion from the core of a young 
stellar cluster. The aims of this study are to determine 1) whether a mass-segregated 
core survives the gas expulsion and 2) the probable location of any massive stars that 
have escaped from the core. Feedback from massive stars is expected to remove the gas 
from the core of the cluster first, as that is where most massive stars are located. We 
find that gas expulsion has little effect on the core for a core star formation efficiency, 
e, of greater than 50%. For lower values of e down to 20%, a reduced core survives 
containing the majority of the massive stars while some of them are dispersed into the 
rest of the cluster. In fact we find that ejected stars migrate from radial to tangential 
orbits due to stellar encounters once they leave the core. Thus, the location of massive 
stars outside of the core does not exclude their forming in the dense cluster core. Few 
massive stars are expected to remain in the core for e lower than 20%. 

Key words: stars: formation - stars: massive - stars: luminosity function, mass 
function - globular clusters and associations: general. 



1 INTRODUCTION 



• The origin of massive stars is a poorly understood phe- 
' nomena. One clue to their formation is that they are gen- 
, erally fou nd in t he dense cores of young st ellar clusters 
jHillenbrand l997i:ICarpenter fc Sanders 1998.: iClarke et alJ 
l200d: iMermilliod fc Garci'al bOOlt IZinnecker et al.1 Il993l) . 
Their location in these very young systems cannot 
generally be explained by dynamical mass segregation 
jBonnell fc D avies 1998). Higher than observed primordial 
core densities (.Kroupa,2002. : Kroupa et al. 200L1 can cause 
rapid mass segregation in the core region. The central lo- 
cation of massive therefore must tell us something about 
their birth. One possibility is that the accretion rates in the 
centre of the cluster are h i gher due to the deepe r potential 
wells (iBonnell et alJl200ll: iBonnell fc Batell2002h and that 
therefore the protostars that are present in the core are able 
to accrete up to a high mass ( eo. lBehrend fc MaedeJ l200lL 
^lessen 2001). An additional possibility is that as the stars 
in the core accrete, the core contracts due to the increase 
in the potential energy. The increase in stellar density may 
then be sufficient such that collisions between intermediate 
mass stars occur to form the massive stars (Bonnell, Bate, 
& Zinnecker 1998; Bonnell fc Ba te 2002). One of the suc- 
cesses of these models is that they naturally account for the 
formation of massive stars in the centres of stellar clusters. 
These models predict the location of massive stars at the 
end of the accretion, mass-buildup phase. At this point the 



cluster contains significant mass in gas which can still affect 
the cluster dynamics. 

Once a massive star has formed in a cluster, it can have 
dramatic efi'ects on its surroundings. The strong winds and 
ionising photons inject significant energy into the surround- 
ing gas which will begin to expand, eventually being cleared 
from the cluster. This gas removal has two effects: firstly, 
it halts any residual accretion or star formation; secondly, 
if the gas still comprises a significant fraction of the total 
mass, the cluster will expand to rcficct the new, lower po- 
tential. Simulations of this process ( Gcycr & Burkcrt 20o3; 
iKroupa et all 1200 it iBoilv fc Kroupal 2003) have concen- 
trated on the global evolution of the cluster as the gas is 
removed. They find that if the gas is a small fraction of the 
total mass, or if it is removed slowly relative to the dynam- 
ical timescale, then the cluster can retain s ignificant frac- 
tion of its initial stars. iKroupa et al.l i200ll) suggests that 
early dynamical mass segregation of the most mass ive stars 
can su rvive during and after a gas expulsion phase. iKroupal 
(^2000) has presented results specifically for models of the 
Orion Nebula Cluster (ONC) including the gas expulsion 
scenario reflecting the situation we model here, but without 
primordial mass segregation. 

In this paper, we are concerned with the dynamical evo- 
lution of the cluster core when the gas is removed. Gas re- 
moval occurs from inside out as the massive stars inject en- 
ergy into their environment. The core of the cluster will be 
cleared first and the core stars will be affected by the gas 
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removal before the rest of the cluster. This has serious impli- 
cations for the dynamics of the massive stars, their location 
in the cluster, or even their possible ejection from the clus- 
ter. We address these questions using numerical simulations 
presented below. 



2 GAS EXPULSION FROM CORE 

Gas removal occurs due to the injection of energy from the 
massive stars in the core. The energy produced by these O 
and B type stars can take the form of strong s tellar winds 
of up to a few x 1000 km s"^ llChurchwelllll999l) or the for- 
mation of an expanding HII region [vhii ~ 10 km s~^). 
In both cases, the expansion velocity is significantly higher 
than the velocity dispersion seen in young stellar clusters 
due to the higher initial density. If the core of the cluster is 
about O.lpc this means that the crossing time for the core 
(tcross ~ 5 X lO^yr) is at least an order of magnitude larger 
than the gas expulsion time. We may then assume that the 
period of gas expulsion has negligible effect on the dynamics 
of the constituent stars in the core, and so we assume instan- 
taneous gas expulsion. Thus we incorporate an initial surplus 
of kinetic energy into the core stars which results in their be- 
ing no longer in virial equilibrium. This implies that the gas 
is only removed from the core as stars further out remain 
in virial equilibrium. The timescales (~ lO^yr) and radial 
extent (~ O.lpc) reflect the situation in an ultra-compact 
HII region where a stellar wind has evacuate d a small in- 
ner region before stalling due to mass l oading (|Lizano_et_aLj 
139t ), ram pressure of the infalling gas (Wood & ChurchwcU 
19891) . or other possible mechanisms ( ChurchwelLi.l999.'l . 
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Figure 1. The mass distribution of the stars in the core. This is 
an average of all initial realizations. This demonstrates the initial 
mass segregation, compared to the non-mass segregated initial 
a = 2.3 distribution shown here as a dotted line. 
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Table 1. Values of the initial gas ratio parameter used, and the 
equivalent star formation efficiency, e. 



3 NUMERICAL SIMULATIONS 

Over the timescales involved it is only gravitational forces 
which play a significant role in the subsequent evolution. 
Thus the system is sufficiently modelled with a pure N-body 
simulation. The numerical simulation s reported her e were 
performed using the N BODY 2 code J AarsethI lioOlh . This 
code is extremely efficient and accurate at following a stellar 
dynamical system for tens to a hun dred dynamical times. It 
includes the AC neighbour scheme iAhmad fc Cohenlll973l) 
and a softened potential (here rgoft ^ 0.1 x rcoro) which min- 
imises computational expense at the cost of not following 
any (close) binary systems. As we are following the stellar 
dynamics over only a few local dynamical times where dis- 
tant interactions dominate, NB0DY2 is sufficient for our 
needs. 

Our cluster is assumed to be in a virialised state prior 
to gas expulsion, from arguments of the star formation 
timescale compared to the crossing times, typical systems 
may indeed be near to viral equilibrium l)Boilv fc Kroup.3 
1200311 . The stellar system consists of 1500 stars populated 
with masses distributed according to a Salpeter mass func- 
tion, 0(m) oc m~^'^, ranging from 0.1 — 10 M© and dis- 
tributed in phase space with a plummer density profile. Ad- 
ditionally, it is necessary to model an initially mass segre- 
gated system such that the massive stars are found in the 
core region of the cluster, as described above. We construct 
such initial conditions by allowing our randomly distributed 



stellar system to relax for several dynamical times until the 
massive stars have sunk to the core. This system is then used 
as our initial condition for the simulations in this investiga- 
tion. It must be emphasised that this process of achieving 
initial conditions is only to construct a model of the ob- 
served initial state and not to represent the physical method 
by which massive stars end up in the cores of very young 
stellar clusters in the first place. Our pre-initial clusters are 
evolved until little more mass segregation can occur. The 
resulting initial mass distribution in the core is shown in 
Figure showing that the high mass stars have sunk to the 
core. 

Our results are presented scaled to the physical units 
which apply to the Orion neb ula cluster as taken from 
iHillenbrand fc HartmannI lll998fl . it should be emphasised 
that our models are not intended to model the ONC which 
has a much wider mass range (0.1 — 50 Mq). The half mass 
radius of the cluster, rhait = O.Spc. We choose the core 
radius to be r-corc = 0.08 pc being 0.1 the half mass ra- 
dius, it is this region within rcoro in which we consider 
the gas expulsion to be effective. The mass of our cluster 
is 1200 Mq, and N= 1500 so that the mean stellar mass, 
< M« >= 0.8 M©. Finally, we follow the evolution of the 
core stars to t > 1.0 Myr 20tcross), by which time the 
dynamical evolution we are interested in will have occurred. 
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Figure 2. Cumulative number ratio for different mass ranges of stars, and different initial gas ratios averaged over all realisations. 
The vertical axis shows the cumulative numbers of stars divided by the total number contained within the region of interest. Radius is 
shown along the horizontal axis. Three mass ranges are shown, indicated by a thick line (M* ^ 5.4 M©), a medium line (5.4 Mq > M« > 
1.08 Mq), and a thin line (M» ^ 1.08 Mq). The dashed lines represent the initial state of the clusters for each mass range. The lower panel 
in each figure shows the cumulative number ratio normalised to the initial numbers of stars within i?core> and the upper panel within 
2-Rcore- The upper figures, (a) Mgas/Mstars = 0.25 and (b) Mgas/Mstara = 1-0, show little evolution; figure (c) Mgas/Matars = 3.0 
shows significant numbers of stars lost from the core; and figure (d) Mgas/^stars = 5.0 show 70% of the massive stars are lost from 
within -Rcorc- 



4 RUN PARAMETERS 

The effect of gas removal is included by way of a super- 
virial distribution of the stellar velocities inside the cluster's 
core. The kinetic energy of these stars is taken to be that 
required to balance an additional mass of gas, Mgas- We have 
parametrised our model by the ratio of Mgas/Mstars prior to 
gas expulsion in the core. The values of this parameter used 



and the equivalent star formation efficiencies are shown in 
Table Q This method mimics the case when gas removal is 
instantaneous. This is the worst case scenario for retaining 
the stars and thus gives a conservative estimate of how many 
massive stars will be retained in the core. 

In practice we take our initial system as described above 
and increase the kinetic energies of the stars in the core by 
the following energy argument. 
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star Formation Efficiency in core 

50% 33% 25% 20% 17% 
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After 0.6Myr, or -10 T^^^^^ (core) 

, I , I 

2 4 

Initial Mg^yM^j^^ in core 



Figure 3. The figure shows the average mass of the 100 stars 
nearest the centre of mass of the cluster for each value of 
Mgas/Mstars- The vertical bars show the resulting range within 
the 8 different realisations of each run. A roughly linear relation 

between Mgas/Mstars and A/lOO after Mgas/Mstars = 1.0. 



[Mgas + Mstars]v1ld = MstarsVn^w 

we can then directly scale the core velocities such that, 



(1) 



(2) 



where we have introduced the star formation efficiency, in 
terms of the gas ratio as. 



1/(1 



Mgas 



(3) 



Mstars 

For each value of the initial parameter, Mgas/Mstars, 8 
initial realisations were constructed using a different random 
seed for each. This helps to extract the mean behaviour of 
the massive stars in the clusters and not individual fluctua- 
tions due to their small numbers. 



5 RESULTS 

The purpose of our investigation is to discover whether a 
mass segregated core will survive after a gas expulsion phase, 
and also to gain general insight to the subsequent motions of 
the massive stars. The following results examine properties 
of the core: mass segregation, average mass within the core, 
fraction of stars remaining in the core, and core survivability. 
We also look at the phase space evolution of the individual 
massive stars by examining their radial evolution and their 
angular momentum. 



5.1 Mass Segregation 

We investigate the evolution of mass segregation in the core, 
which, in this case, amounts to the survival of the most 
massive stars in the core. 

After 0.62 Myr (~ 12fcross) we take the ratio of the num- 
ber of stars within r to the number initially contained within 
given radii (i?corc and 2i?corc). This cumulative number ra- 
tio is calculated independently for 3 stellar mass ranges: 



M. 1.08 Mq, 1.08 M0 < < 5.4 M0, 5.4 M© ^ M*, 
shown as a thin line, a medium line, and a thick line respec- 
tively in Figures 121 a)-(d). These figures show the cumula- 
tive number ratio for a Mgas/Mstars parameter of (a) 0.25, 
(b) 1.0, (c) 3.0, and (d) 5.0. The top and bottom panels in 
each figure show the cumulative number ratio normalised to 
27?core and Rcore respectively. The initial state of the sys- 
tem at f = is shown as a dashed line in each figure. These 
figures are the combined results of all realisations for each 

particular Mgas/Mstars- 

In particular, we are interested in the highest mass 
range, which we consider to represent the OB type stars 
in the cluster. Clusters with lower values of Mgas/Mstars 
(figures I^I a) and (b)) show little significant change in the 
distribution of the massive stars, in fact the lower mass stars 
decrease in number more than the massive stars, making the 
core slightly more segregated. However, as soon as the effec- 
tive kick to the core stars becomes large enough to boot out 
the massive stars we see a dramatic decrease in all the stellar 
numbers. Figures |5Jc) and (d) show a decrease in the actual 
massive star numbers by 20% and 70% respectively. 

An alternative and perhaps more intuitive way of look- 
ing at the mass segregation is to look at the average stel- 
lar mass within a certain radius, or as we do in this case 
the 100 stars nearest to the centre of mass. Figure |3 shows 
there is an almost linear decrease in Mioo for each value of 
the Mgas/Mstars parameter after Mgas/Mstars = 1-0. Also 
shown, as a vertical bar on each point, is the range due to the 
variation in the initial realisations for each parameter. This 
shows clearly that as the stars are ejected out into the rest 
of the cluster, the massive stars are less of a characteristic of 
the central cluster regions. Although the mean stellar mass 
is slightly greater than the cluster mean < A/, 0.8 M©, 
at Mgas/Mstars = 5.0 the mass segregation is much less 
significant. 

5.2 Radial evolution of Massive stars 

The question now arises of what becomes of the massive 
stars which the core has ejected. For comparison with ob- 
served clusters we would like an idea of where the stars are 
located in our simulations. We examine the radial distances 
of the most massive stars from the centre of mass with re- 
spect to time. Typical examples are shown in Figure 2] 

These figures illustrate typical evolution of the 10 most 
massive stars which were initially in the cluster core, each 
figure showing one realisation for each value of the parameter 
Mgas/Mstars- As wc have seen in the previous section, for a 
larger Mgas/Mstars, & greater number of massive stars leave 
the core. Figure |l]shows this trend explicitly and we can see 
the subsequent radial motion of the stars after gas expulsion. 
The grey region in the figures is the initial core in which the 
gas is assumed to have been removed. The half mass radius, 
i?haif, is also shown as a solid line at 0.8pc. 



For a Mgas/Ms, 



0.25, the figure shown top left is 



typical for all realisations and similar results are seen for 

Mgas/Mstars = 1.0 (tOp middle). For Mgas/Mstars = 2.0 

(top right) it is more typical to see a couple of the stars 
ejected from the core but within the half mass radius. For 
Mgas/Mstars ~ 3.0, 2 to 4 stars are ejected beyond JJhaif 
and another couple beyond Rcoic (bottom left). Figure |^ 
shows a resulting configuration in the x — y plane. There 
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Figure 4. Radial distance from density centre vs time for the 10 most massive stars initially found in the cluster core. These figures 
represent a single individual realisation for each value of the parameter Mgas/Mstars- The grey area is the core region in which the 
initial velocity kick is given, and a solid line representing the initial -Rhalf is shown. The thickness of the individual lines relate to the 
mass of the star, thickest being most massive. The realisations have been chosen to illustrate the typical behaviour of the system for 



each M„ 



where no stars leave the core for a value of 0.25 and all leave the core for Mgaa/Ms. 



■ 5.0. The apparent kink in 



the path of stars that leave the core is due to the time sampling of the simulations and is not a physical affect. The value of the lower 
mass limit of stars shown is also indicated in each figure. 



0.62Myr 




X — pc 



Figure 5. The x — y plot of a typical realisation with 
Mgas/Mstars = 3 after t = 0.62 Myr (~ 12icross). Here we notice 
two stars remain bound within -Rhalf (large dashed circle) and the 
rest forming a bound (but less dense) core shown in the inset. 



is a resemblance to the configuration of the massive stars 
in the Orion Nebula Cluster (ONC), showing that this is a 
possible mechanism leading to ONC type structures, where 
some of the massive stars are located a significant distance 
from the core. Although without observational velocities we 
cannot say how likely it is. 

The tendency when Mgas/Mstars ~ 4.0 (bottom mid- 
dle) is for more massive stars to be completely ejected 
and still two or so to be found between -Rcore and /?haif. 
It is likely tliat most stars are ejected from the core for 
Mgas/Mstars ~ 5.0 and also beyond i?haif (Figure 0| bot- 
tom right) . 

5.3 Velocity Evolution 

Observational measurements of velocities (proper motions as 
well as line-of-sight velocities) are now feasible with space 
borne instruments, such velocities may provide strong evi- 
dence to decide on the spatial origin of the stars. It is there- 
fore of interest to explore the kinematics that result from gas 
expulsion and the ejection of massive stars from the core of 
the cluster. 

The expectation is that stars which receive significant 
velocity increases escape the cluster on near-radial orbits. In- 
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stars 



almost exclusively increase their angular momenta such that 
1000 > |Lt|/|-Lo1 > 1. 

Such a result indicates that the higher velocity massive 
stars must have encountered a greater number of stars with 
a broader range of orbits on their journey outward through 
the cluster. A combination of their gain in \L\ and high 
velocity puts many stars on large tangential orbits which do 
not pass through the core. This was seen in Figure^Jas stars 
are ejected from the core but are seen to orbit beyond Rcoib- 
As a consequence of this the velocities of the ejected massive 
stars will be seen to have predominantly tangential motions 
to the cluster core. 

There is a further subset of stars which appear to leave 
the cluster entirely (a limit of 5i?haif is used to identify 
cluster leavers), these are shown as open triangles in Fig- 
ure |S| The physical velocities of these stars range from 10 
to 40kms~^. A result which allows us to rule out the ex- 
istence of observed high velocity stars 200kms~^) as a 
consequence of the two body interactions which occur in this 
scenario. 



Figure 6. Ratio of angular momentum \Lt\ at i = 
0.62 Myr(12tcross) to initial angular momentum \Lo\ for all stars 
over 5 Mq. All the data for each realisation of each M gas /M stars 
is shown, separated into stars which are still in the core at 
t = 0.62 Myr (o, shifted left for clarity), stars beyond the core 
(*), and stars which have left the cluster (A shifted right, for 
r > 5-Rhaif)- The distribution shows a distinct bimodality be- 
tween core stars which show little overall increase in \L\, and 
stars that have left the core which all show an increase in \L\ 
becoming greater with higher M gas /I^I stars 



stead, we show below that these stars generally have more 
tangential than radial motions. This can be explained by 
the interactions that occur due to the discrete nature of the 
stellar cluster. We can suppose that the orbits of these stars 
will encounter a great range of stars on highly varied orbits, 
its velocity signature will vary greatly, leaving the angu- 
lar momentum measured with respect to the centre of the 
core density peak, \L\, as the most robust general measure 
of the stars' subsequent evolution. This measure of angular 
momentum is not affected by the motion of core over the 
duration of these simulations. 

Figure |S| shows the ratio of the angular momentum at 
time t — 0.62 Myr(12fcross) to the initial value immediately 
after the velocity kick \Lt\/\Lo\. All stars over SM©, in all 
realisations of the parameter Mgas/Mstars are shown. We 
have distinguished the stars by their radial position at time 
t — 0.62 Myr: those within -Rcore are shown as small open 
circles (o), those beyond -Rcore shown as filled stars (*), and 
those beyond 5-Rhaif (those which have left the cluster) are 
shown as open triangles. Each radial range is shifted slightly 
to left or right for clarity. 

The first effect to notice is the increasing spread of an- 
gular momentum with increasing initial Mgas/Mstars, such 



that for Mgas/Ms, 



5.0 the gain in \L\ for some stars 



reaches ~ 1000. The other striking feature is the bimodal- 
ity between stars that have stayed in the core and stars 
which have left it. Stars remaining in the core have angular 
momenta distributed around their initial values, although 
plenty of angular momentum exchange has occurred result- 
ing in 100 > lirl/liol > 0.01. The core escapers, however. 



5.4 Core Survival 

To summarise the results in terms of the original ques- 
tion posed — how does gas expulsion from the core affect 
the core's survival? — we have plotted (Fig the frac- 
tion of massive stars remaining in the core for each value 
of Mgas/Mstars- It is pcrhaps a matter of definition as to 
what constitutes a surviving core, and in this case we take a 
value of 50% of the original massive stars. More precisely we 
say that if more than half the core (or rather its constituent 
massive stars) is dissipated subsequent to the gas expulsion 
phase then the core does not survive. This is supported by 
the fact that the actual number of massive stars in the core 
are of the order of a few, once half of these have left the 
system there may be only one or two left, not enough to 
constitute a core. 

The relation between Mgas/Mstars and the core sur- 
vivability is fairly smooth (figure [7|. The core is relatively 
undisturbed up to Mgas/Mstars ~ 2.0, then becomes pro- 
portionately more disrupted with greater Mgas/Mstars but 
survives (by our definition) to a Mgas/Mstars of over 4.0. 



6 DISCUSSION 

We have modelled the gas expulsion from the core of a stellar 
cluster via a rapid change in the stellar velocities and thus 
an effective change in the gravitational potential of the core. 
This is taken to mimic the evolutionary stages of a hyper 
or ultra-compact HIT region. Their evolution is likely to be 
driven by a stellar wind which stalls as mass loading occurs 
in order to explain their relatively long lifetimes (~ lO^yr). 
Thus, we assume an instantaneous removal of gas from the 
core of the cluster but assume that the gas content of the 
rest of the cluster is for the most part unaffected during our 
simulation. 

The formation of massive stars is largely an unsolved 
problem. What is certain is that they require high accretion 
rates in order to build the star up over timescales of 10^ to 
10^ yea rs. This is true whether they form just through ac- 
cretion iBehrend fc MaedeillioOll : iBonneU fc Batell2002ll or 
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star Formation Efficiency in core 

80% 50% 33% 25% 20% 17% 



interactions with a binary or multiple system which eject one 
or more of the stars. Some of these ejections could be the 
high-velocity OB runaways while those that do not garner 
such a large velocity boost remain in the cluster but outside 
the core. If that is the case, their dynamics should resemble 
the massive stars ejected from the core in the simulations 
reported here. 



After 0.6Myr, or ~10 T (core) 

^ cross^ ^ 



Figure 7. Quantifying core dissolution for different initial 
Mgas/Mstars values, and the equivalent star formation efficien- 
cies, after a period of 0.62Myr {12tcross) 



even through stellar mergers which require high accretion 
rates to reach the necessary stellar density (j^onncU ct al. 
Il998l: iBonnell fc BatdbOOsli . The most obvious way that 
high accretion rates can occur is by situating the forming 
massive star in the centre of a stellar cluster. In this way, 
the overall cluster potential funnels matter down to the pro- 
tomassive star where i t can preferentially accrete it due to 
its position and mass llBonneU et alJl2001^ . This naturally 
results in a mass segregated stellar cluster even while deeply 
embedded. 

The simulations reported here show that gas expulsion 
is unlikely to remove this 'initial' mass segregation unless the 
gas comprises more than 83 per cent of the total mass in the 
core. The core, being where the accretion rates are highest 
and the most massive stars are located, has therefore the 
lowest gas fraction of the system and is the most likely to be 
dominated by the mass in stars. Thus, the mass segregation 
which results from the cluster formation and gas accretion 
is likely to remain unaffected by the gas expulsion stage in 
the cluster. 

In some young clusters such as the ONC, a few rela- 
tively massive stars are found to be some distance from the 
core of the cluster. If the above scenarios for the forma- 
tion of massive stars are correct in that they form in the 
denser parts of the cluster, then we need to explain how 
they could have moved from the core to, for example, the 
half-mass radius of the cluster, the numerical simulations 
reported here offer one possibility. The massive stars could 
have originated in the core but escaped from this region dur- 
ing the gas expulsion phase (Williams 2002). We see that for 
gas fractions of greater than 83 per cent, some of the massive 
stars do escape the core and can be found virtually anywhere 
in the cluster. They will of course sink back to the cluster 
core but on longer timescales than the age of the cluster 
(BonncU fc Davios 199^. 

The difficulty in the above is that relatively high gas 
fractions are required which may not be present in the core 
of the cluster. An alternative mechanism to explain the lo- 
cation of the massive stars not in the core would be stellar 



7 CONCLUSION 

Using a series of N-body simulations we follow the dynami- 
cal evolution of massive stars which are assumed to be ini- 
tially located in the core of the cluster. The stars undergo 
evolution because of very rapid gas expulsion (assumed in- 
stantaneous) which is parametrised by the ratio of gas mass 
to stellar mass in the core (Mgaa/Mstars)- This parameter 
is varied from 0.25 to 5.0 (equivalent to a star formation 
efficiency of 80% down to 17%). 

We find that the initial mass segregation in the core is 
destroyed for an initial gas mass of greater than 4 times the 
stellar mass, or a star formation efficiency of less than 20%. 
The initial conditions used are of a relaxed, mass-segregated 
cluster. The massive are therefore approximately in equipar- 
tition with low-mass stars in the core. Their lower velocities 
result in a lower probability of escaping the core once the 
gas too has been removed. This explains the relatively high 
value of Mgas/Mstars required to dissolve the core. Other- 
wise, for star formation efficiencies of more than 20%, but 
less than 50% some massive stars will be ejected from the 
core and maybe beyond. However these ejected stars are 
found to have predominantly tangential motions, gaining 
angular momentum due to encounters along the way. 
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